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The magnetic interactions in one-dimensional, two-dimensional (2D) and ladder cuprates are evaluated system- 
atically by using small Cu-0 clusters. We find that the superexchange interaction J between nearest neighbor 
Cu spins strongly depends on Cu-O structure through the Madelung potential, and in 2D and ladder cuprates 
there is a four-spin interaction J cyc , with magnitude of 10 % of J. We show that J cyc has a strong influence on 
the magnetic excitation in the high-energy region of 2D cuprates. 



A variety of insulating cuprates affords us an 
opportunity to study the magnetic properties of 
low-dimensional systems. Recent experiments 
for insulating cuprates have revealed interesting 
characteristics of magnetic interactions: The su- 
perexchange interaction between nearest neigh- 
bor Cu spins J remarkably depends on Cu-0 
network structure JjJ , and additional interactions 
such as a four-spin (4S) interaction are important 
for ladder |2| and two-dimensional (2D) cuprates 
|||. These characteristics indicate the neces- 
sity to establish proper magnetic descriptions for 
the cuprates. In this paper we perform a sys- 
tematic study of magnetic interactions for one- 
dimensional (ID), 2D, and ladder cuprates theo- 
retically. 

A starting model to describe the electronic 
states of cuprates is the d-p model, in which hop- 
ping integrals between Gu3d and 02p orbitals 
(T p d) and between 02p orbitals (T pp ), an energy- 
level separation between the Cu3rf and 02p or- 
bitals (A), and Coulomb interactions on Cu and 
O sites are taken into account. T p d and T pp are 
obtained by considering not only the bond length 
dependence but also the effect of the Madelung 
potential around Cu and O ions. We find that 
the potential enhances the magnitudes of T pc i and 
T pp in the ID cuprates as compared with those in 
the 2D ones . In the two- leg ladder compounds 
such as SrCu203, T pp along the leg of the ladder 
is enhanced by the Madelung potential due to ad- 
jacent two-leg ladders. These enhancements play 
an important role in the dependence of J on the 



dimensionality. The A is determined from the 
difference in the Madelung potential between Cu 
and O sites. 

The magnetic interactions are evaluated by 
mapping the lowest several eigenstates of small 
clusters (CU2O7, CU4O12, and CU6O17) for 
the d-p model onto those of the corresponding 
Heisenberg-type model For 2D systems, we 
take into account not only J, but also a diagonal 
interaction Jdi ag and 4S interaction J cvc in the 



model: H = J£ (i ,-vS 
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p ijki)> where s * is a 



i>j + Jcyc Zjplaquettei "ijkl + 

spin operator at site i, and J cyc is defined as the 
coefficient of the 4S cyclic permutation operators 
Pijki and Pj~M , which can be rewritten by using 
the two-spin interaction (Sj - Sj) and the four-spin 
interactions (Si • Sj)(Sk - Si). For ladder systems, 
we distinguish between the nearest neighbor in- 
teractions along the leg ( Ji cg ) and along the rung 
( Jrung) of the ladder. 

The calculated results are summarized in Ta- 
ble 1, where we take La2Cu04, SrCu203, and 
Sr2Cu03 as typical systems of 2D, ladder and 
ID cuprates, respectively (see Refs. 1 and 5 for 
the parameters used in the calculations). We find 
that J in the ID cuprate is larger than that in 
the 2D one. This is caused by the enhancement 
of the hopping integrals in ID cuprates as men- 
tioned above. For 2D cuprates, we obtain J to be 
~0.15 eV, consistent with the experimental val- 
ues 0. In addition, we find that J cyc is 7 % of J, 
while Jdiag is zero. These results are consistent 
with a previous cluster calculation [pi, and a re- 
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Table 1 

The magnetic interactions obtained by mapping 
the eigenstates of CU2O7, CU4O12, and CU6O17 
clusters for Sr 2 CuC>3, La2Cu04, and SrCu203, 
respectively, onto the corresponding Heisenberg- 
type models. The numbers in parentheses denote 
the deviation in the last significant digit. 



Material 


J[cV] 


Jdiag [eV] 


Jcyc [cV] 


Sr 2 Cu0 3 


0.17 






La 2 Cu04 


0.146(1) 


0.00(0) 


0.011(1) 


SrCu 2 3 


Ji og : 0.195(5) 


0.003(2) 


0.018(2) 




J rung : 0.15(2) 
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Figure 1. The magnetic excitation spectra of a 
4x4 Heisenberg model with J and J cyc . 



cent analysis of a multimagnon spectrum J}] . For 
ladder cuprates, we obtain Ji eg /Jrung = l-3- The 
enhancement of T pp along the leg of the ladder is 
the origin of the relation Ji eg > Jrung- Jcyc is 10 % 
of Jieg. Note that the ratio Ji eg / Jrung is smaller 
than that believed so far, i.e., Ji e g/Jrumg~2 (J. 
However, the Heisenberg ladder with the present 
values of Ji eg , J rU ng, Jcyc and Jdi ag reproduces 
very well the experimental results of the temper- 
ature dependence of the magnetic susceptibility 
(not shown here). Therefore, we consider the val- 
ues shown in Table 1 to be reasonable. Here we 
would like to emphasize that J cyc plays a crucial 
role in obtaining the good agreement between the 
experimental and theoretical magnetic suscepti- 
bilities J§. 

Next, in order to examine the effect of J cyc 
on the magnetic excitation, we calculate the dy- 
namical spin-correlation function S(q,u>) for 2D 
cuprates. Figure 1 shows the dispersion and the 
intensity of S(q, w) for a 4x4 Heisenberg model 
with J=0.146 eV and J cyc =0.011 eV. For compar- 
ison, the result for J cyc =0 are also shown. We find 
that the intensity is not sensitive to J cyc , while 
the dispersion is strongly suppressed in the high- 
energy region. In particular, it is worth noting 
that w(q) at q=(7r/2,7r/2) becomes smaller than 
that at q=(7T,0). This is in contrast with the case 
of J cyc =0, in which the magnetic zone bound- 
ary (q=(7r,0)— >(7r/2,7r/2)) has a flat dispersion. 
Therefore, it is desirable that inelastic neutron- 
scattering experiments in the wide energy region 



be performed to verify the role of J cyc , that is, 
the suppression of the dispersion at (ir/2,n/2). 

In summary, we have evaluated the magnetic 
interactions in various cuprates systematically. 
We have shown that an ionic nature inherent in 
insulating cuprates is important for the material 
dependence of J. We found that J cyc is ~10 % of 
J, and greatly influences the magnetic excitation 
spectra in 2D cuprates. 
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